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ABSTRACT 

The  transient  excitation  of  a cavity  through  an  aperture  is  used  to 
illustrate  the  state-space  solution  of  transient  electromagnetic  problems. 
Depending  on  whether  the  E-  or  the  H-formulation  is  used,  the  effect  of 
the  aperture  can  be  accounted  for  by  a magnetic  current  in  invoking  the 
induction  theorem  or  by  an  electric  current  in  invoking  the  equivalence 
theorem.  The  governing  second-order  differential  equation  is  converted 
into  a set  of  first-order  state  equations  by  defining  three  new  state 
variables  in  addition  to  an  appropriate  vector  potential.  These  state 
equations  are  solved  by  the  method  of  moments.  Two  cases  are  considered: 
parallel  polarization  and  perpendicular  polarization.  Some  significant 
singularities  for  the  parallel-polarization  case  are  found  and  the  elec- 


tric intensities  as  functions  of  time  at  two  locations  in  the  cavity  are 


computed  for  a step  excitation 


The  Findings  in  this  report  are  not  to  be  construed  as  an  official  n, 
Department  of  the  Army  position,  unless  so  designated  by  other  author 
ized  documents. 
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I . INTRODUCTION 


Transient  electromagnetic  problems  arise  in  the  consideration  of 
shielding  effectiveness,  electromagnetic  compatibility,  radar  target 
identification,  remote  environmental  sensing,  and  electromagnetic  effects 
due  to  X rays  and  gamma  rays  generated  by  a nuclear  blast.  The  time- 
harmonic  behavior  of  an  antenna  or  a scatterer  can  be  obtained  first  as  a 
function  of  frequency  by  solving  the  governing  integral  equations  with  the 
given  boundary  conditions.  The  transient  response  can  then  be  determined 
by  performing  an  inverse  Fourier  transformation.  Numerical  methods  are 
used  in  the  frequency-domain  solution  for  many  different  frequencies.  The 
brute-force  numerical  Fourier  inversion  is  generally  inefficient  and  con- 
vergence problems  arise  in  superposing  the  steady-state  solutions.  A few 
simple  situations,  such  as  transient  scattering  from  wire  antennas  and  con- 
ducting cylinders  [l]-[3],  have  been  analyzed  directly  in  the  time  domain. 
Basis-function  expansions  and  inner  products  over  both  space  and  time  are 
required.  Space-time  integro-dif ferential  equations  are  encountered  and 
the  numerical  representations  of  the  derivatives  lead  to  very  complicated 
procedures . 

A variation  of  the  time-domain  approach  makes  use  of  a a Hall^n-type 
integral-equation  formulation  [4].  As  a consequence  of  the  absence  of  space- 
time  derivatives  under  the  integrals,  the  numerical  process  is  less  compli- 
cated. However,  because  of  retarded  time,  a tedious  step-by-step  iterative 
procedure  has  to  be  used  in  conjunction  with  the  proper  families  of  trajec- 
tories known  as  characteristic  curves  in  order  to  determine  the  homogeneous 
solutions.  All  of  the  methods  mentioned  above  are  tedious  and  do  not  afford 


a physical  insight  in  the  solution  of  a transient  problem.  Moreover,  a 
complete  recalculation  would  be  necessary  under  any  change  in  the  wave 
shape,  polarization,  or  angle  of  incidence  of  the  source  of  excitation. 

More  recently  a singularity-expansion  method  has  been  used  to  deter- 
mine the  transient  response  of  scatterers  of  simple  geometries  [5 ] — f 11 ] . 
Responses  to  transient  excitations  are  expressed  in  terms  of  exterior 
natural  frequencies,  modes,  and  coupling  coefficients,  and  induced  currents 
are  represented  by  a series  of  damped  sinusoidal  functions.  This  method  has 
the  advantages  of  providing  a physical  sight  to  the  radiation  or  scattering 
problem  and  of  allowing  the  response  to  be  determined  for  a change  in  source 
parameters  without  a complete  recalculation.  In  the  evaluation  of  the  natural- 
mode and  coupling  vectors,  it  is  necessary  to  know  the  nature  of  the  singu- 
larities involved  in  the  inverse  Laplace  transform.  Knowledge  in  this  respect 
is  not  yet  secure.  The  transcendental  nature  of  its  system  impedance  matrix 
results  in  an  infinite  number  of  complex  poles  whose  locations  must  be  numeri- 
cally searched.  It  is  generally  agreed  that  the  singularity-expansion  method 
is  not  very  satisfactory  for  evaluating  early-time  responses. 

The  method  of  characteristic  modes  has  been  used  for  determining  the 
steady-state  response  of  conducting  bodies  [12].  A striking  similarity  appears 
to  exist  between  this  method  and  the  singularity-expansion  formulation  if  the 
method  of  moments  [13]  using  a common  spatial  basis  is  applied  to  both  cases; 
but  no  formal  relations  have  yet  been  established.  The  natural  modes  from 
the  singularity  expansion  are  not  orthogonal.  The  characteristic  modes  are 
orthogonal,  but  they  vary  with  the  source  frequency  and  it  would  be  necessary 
to  compute  the  characteristic  modes  for  all  frequencies  before  a Fourier  inver- 
sion could  be  effected  to  determine  the  transient  behavior. 


The  problem  of  transient  field  behavior  inside  a conducting  cavity 
due  to  excitation  through  an  aperture  by  an  incident  electromagnetic  pulse 
(EMP)  is  particularly  difficult  because  of  the  reflections  from  the  cavity 
walls  and  the  coupling  between  the  interior  and  exterior  fields  at  the  aper- 
ture. Past  investigations  on  EMP  excitation  of  cavity-backed  apertures  have 
largely  dealt  with  small  openings  and  have  neglected  the  effect  of  cavity 
reflections  on  the  aperture  field  distribution.  For  small  openings  the 
quasi-static  method  is  used  to  determine  the  fictitious  magnetic  current 
and  charge  distributions  in  the  aperture.  Equivalent  electric  and  magnetic 
dipoles  are  defined,  and  their  radiated  fields  determined  with  the  aid  of 
scalar  and  vector  potentials  [ 14 ] — [ 19 ] . The  fields  in  the  cavity  are  custom- 
arily expanded  in  terms  of  unperturbed  normal  modes.  The  quasi-static  approxi- 
mation cannot  be  applied  when  the  aperture  is  not  small  and  when  early-time 
responses  are  important.  In  neglecting  the  effect  of  the  reflections  from 
cavity  walls  on  the  aperture  field  distribution,  one  essentially  treats  the 
external  and  internal  portions  of  the  problem  separately.  Since  cavity  dimen- 
sions obviously  play  an  important  part  in  the  total  problem,  this  approach  may 
result  in  significant  errors. 

In  this  report  we  will  avoid  the  quasi-static  approximation  and  solve 
the  internal  and  external  portions  of  the  problem  simultaneously.  New 
variables  (state  variables)  will  be  introduced  to  convert  the  governing 
second-order  differential  equation  into  a set  of  first-order  equations  which 
correspond  to  normalized  state  equations.  The  field  within  the  cavity  will 
be  expanded  in  terms  of  suitably  chosen  subsectional  expansion  functions  with 
variable  coefficients  and  the  field  outside  the  cavity  expressed  as  a super- 
position of  plane-wave  fields.  TTie  cavity  and  the  external  fields  are  matched 
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at  the  aperture  where  a fictitious  equivalent  current  exists.  A combined 
field  expression  containing  the  unknown  expansion  coefficients  is  obtained. 

To  determine  these  coefficients  the  method  of  moments  [13]  is  used  to  con- 
vert the  first-order  equations  into  matrix  equations.  It  will  be  shown  that 
the  typical  coefficient  matrix  can  be  expressed  in  a form  for  which  the 
singularity-expansion  method  [5]  can  be  used  to  advantage. 

The  general  procedure  of  solution  for  cavity-backed  aperture  problems 
is  outlined  first.  The  theoretical  formulation  for  the  transient  excitation 
of  a rectangular  cavity  with  a slot  aperture  is  then  given  for  both  parallel 
and  perpendicular  polarizations.  Some  numerical  results  are  included. 
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II.  SOLUTION  PROCEDURE 


We  consider  the  problem  of  a slot  aperture  in  an  infinite  conducting 
plane  backed  by  a rectangular  conducting  box,  as  shown  in  Fig.  1.  An  inci- 
dent transient  electromagnetic  wave  (E^ , H^)  impinges  normally  on  the  plane 
and  the  aperture.  The  problem  is  to  determine  the  scattered  field  in  the 
y > 0 region  and  the  field  penetrated  through  the  aperture  into  the  conduct- 
ing cavity. 

The  induction  theorem  [20],  [21]  can  be  invoked  for  the  solution  of 

this  problem.  Figure  2(a)  represents  a simplified  2-dimensional  view  of 

the  original  problem.  (E  , H ) and  (E  , H ) are,  respectively,  the  cavity 

c c s s 

field  and  the  external  scattered  field.  In  order  to  determine  these  unknown 

fields,  we  consider  the  case  when  the  aperture  is  covered  by  a conductor.  The 

entire  region  to  the  left  of  the  infinite  plane  will  have  a null  field  and, 

according  to  the  induction  theorem,  a magnetic  current  on  the  right  surface 

of  the  conducting  plane  will  support  a different  scattered  field  (E°,  H°) , as 

s s 

shown  in  Fig.  2(b),  where 

M = E°  * n 
o s 

= fl  x E1 

= f x E1  . (1) 

For  a normally  incident  plane  wave  (E*,  H*) , the  scattered  field  (Eg,  Hg) 

from  an  infinite  conducting  plane  without  an  aperture  is  easily  determined. 

The  null  field  to  the  left  of  the  plane  will  be  maintained  if  the  plane  is 

removed  and  a magnetic  current  2Mq  exists  in  its  place  which  will  result  in 

a field  (E*  + E°,  H1  + H°)  in  the  y > 0 region,  as  shown  in  Fig.  2(c). 
s s 

I 
i 

| 
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Subtracting  the  fields  in  Fig.  2(c)  from  those  in  Fig.  2(a),  we 
obtain  the  problem  in  Fig.  2(d).  The  magnetic  current  M in  the  aperture 
Is 

M = . 2M  = - 2y  x E1  , (2) 

o 

which  supports  the  field  (E  , H ) inside  the  cavity  and  a field  (E  - E°, 

c c s s 

H - H°)  to  the  right  of  the  infinite  plane.  We  note  that  the  region  in 
s s 

which  the  difference  field  (E  - F.°,  H - H°)  exists  is  source-free  and 

s s s s 

that  the  tangential  component  of  the  electric  field  is  required  to  vanish 
on  conducting  walls. 

For  the  problem  in  Fig.  2(d),  we  start  from  the  two  Maxwell's  curl 
equations 

V x E = - v - M (3) 

O dt 

V x h = e ||  . (4) 

O dt 

Taking  the  Laplace  transform  of  Eqs.  (3)  and  (4),  we  obtain 

V*E=-tJsH-M  (5) 

o 

V x H = c s 1 . (6) 

o ' 

where  a tilde  (-)  over  a quantity  denotes  the  Laplace  transform  of  that 
quantity. 

Let  F be  the  Laplace  transform  of  an  electric  vector  potential  F 
such  that 

E = - V x F . (7) 

Combining  Eqs.  (5)  to  (7)  and  using  the  Lorentz  gauge,  we  have  an  in- 
homogeneous Helmholtz  equation: 
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V2F  = U € S2  F = - M . (8) 

o o 

Solution  of  Eq.  (8)  for  F will  give  E from  Eq.  (7)  and  H from 

H = — V x V x F (9) 

V 

in  regions  where  M is  zero. 

Specialization  of  these  general  formulas  will  depend  on  the  polari- 
zation of  the  incident  wave;  but  as  soon  as  E*  is  known,  Eqs.  (7)  and  (9) 
can  be  expanded  into  component  equations  and  the  source  term  M in  Eq.  (8) 
can  be  found  from  Eq.  (2).  A set  of  new  variables  can  then  be  defined 
which  will  convert  the  second-order  differential  equation  (8)  into  a set  of 
first-order  equations,  and  these  equations  are  Laplace-transformed  normalized 
state  equations  in  the  new  state  variables.  Solution  of  the  transformed 
state  equations  involves  four  steps.  First,  the  space  inside  the  cavity  is 
divided  into  subsections  and  suitable  expansion  functions  are  chosen  over 
the  subsections.  The  elements  of  the  unknown  state-variable  vector  are  then 
expressed  in  terms  of  the  expansion  functions  within  the  cavity.  Second,  the 
field  in  the  y > 0 region  is  expressed  as  a superposition  of  plane  waves. 
Third,  the  cavity  and  the  half-space  fields  are  matched  at  the  aperture. 
Fourth,  inner  products  are  taken  so  that  the  matrix  equations  for  the  un- 
known expansion  coefficients  are  obtained.  These  steps  are  outlined  sepa- 
rately for  the  cases  of  parallel  and  perpendicular  excitations  in  following 
sections . 


w 


In  this  section  we  consider  the  case  of  an  incident  plane  wave  with 
the  electric  field  polarized  in  a direction  parallel  to  the  slot.  Referring 
to  Fig.  1,  we  have 

E1  = 2 E1  (10) 

z 

and  the  Laplace  transform  of  Eq.  (2)  becomes 

M = - 2H  E1  = ft  M (11) 

z x 

which  has  only  an  x-component.  The  x-component  of  Eq.  (8)  is  then 


~2~  2 - 
V F -yes  F 

X o O X 

where  6(y)  is  a Dirac  delta  function. 


= - M 6(y) 
x > 

From  Eqs.  (7)  and  (9), 


(12) 


we  have 
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y s 3x3y 
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y s 3x3z 
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F . 
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(13) 

(1A) 

(15) 

(16) 

(17) 

(18) 


The  second-order  differential  equation  (12)  can  be  represented  as  a 
set  of  first-order  equations  by  defining  new  quantities  u,  v,  and  w such 


(19 


- Fx(r»s)  = s a(r,8) 

- Fx(r,s)  = s v(r,s)  (20 

and 

- 4~  F (r,s)  = s w(r,s)  , (21 

O Z X 

where  r Is  the  space  variable.  We  have,  from  Eq.  (12), 

4 Q(r’s)  + Hr’S)  + h *(r>s)  = ' PoeoS  fx(r»s)  + 7 \ ■ (22 

Comparing  Eqs.  (21)  and  (20)  with  Eqs.  (14)  and  (13)  respectively,  we  see 
that 

E = s w (23 

y 

and 

E = - s v . (24' 

z 

The  introduction  of  u,  v,  and  w and  the  use  of  the  first-order  equations 
will  result  in  significantly  faster  convergence  in  the  numerical  solution 
The  first-order  equations  (19)  to  (22)  can  be  written  in  a succinct 
form  by  defining  the  following  operators  and  column  matrices: 


L f (r,s)  = s P f(r,s)  + e^(r,s) 


(29) 


Note  that  the  Inverse  Laplace  transform  of  Eq.  (29)  Is  a set  of  normalized 
state  equations  in  the  four  state  variables  F , u,  v,  and  w. 

In  order  to  solve  Eq.  (29)  by  the  method  of  moments  we  subdivide 
the  space  within  the  cavity  in  the  x,y,  and  z directions  and  choose  expan- 
sion functions  Fx(1>j  §k)  <r)  , “(i.j.k)00’  v(i,j,k)(r)’  and  W(i, j ,k) (r)  0Ver 
the  subsections.  The  expansion  functions  must  satisfy  the  required  boundary 
conditions.  For  convenience,  we  define  the  following  column  vectors: 
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f(i.J.W<r) 


x(i,j,k)vw 

0 

0 

0 

0 

0 


v<i,j,k><r) 


(i.j.k) 


(r)  = 


£(i,1,w(r)  ■ 


u(i>j  .k) £t) 
0 

0 


0 

0 

0 


v<i.j.k)<rt 


(30) 


In  view  of  Eq.  (27),  we  can  then  write  the  expanded  form  of  f(r,s)  inside 
the  cavity  as 


l(r,s)  = l {a(J  J ,.N(s)  f*4  uN(r)  + (s)  f“,  4 


i.j.k 


‘(i.j.k) l(i,j,k)vw  T p(i,j,k)^  l(i,j,k) 


+ Y(i,j,k)(s)  f(i,.1,k)(r)  + 6(i,j,k)(s)  f(i,j,k)(r)}' 


Note  that  the  expansion  functions  F , u,  v,  and  w are  functions  of  position 
only  and  that  the  inverse  transformation  of  5,  (3.  Y,  and  5 will  yield  the 
time-varying  expansion  coefficients. 

The  field  in  the  half-space  y > 0 is  expressed  as  a superposition  of 
plane  waves  and  the  internal  and  external  fields  are  matched  at  the  aperture 
y * 0.  A combined  field  expression  for  f(r,s)  can  be  obtained  which  holds 
Inside  the  cavity,  in  the  y > 0 half-space,  as  well  as  in  the  slot.  This 
is  then  substituted  in  Eq.  (31)  which,  after  inner  products  with  the  expan- 
sion functions  in  Eq.  (30)  have  been  taken,  leads  to  the  following  matrix 


equation: 
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where  m , . , , are  expansion  coefficients  for  the  magnetic  current  M , 

x v 1 y j x 

C(s)  is  a column  matrix  and  J l 's,  p 's  and  q's  are  themselves  matrices 

um  mn 

arising  from  the  inner  products,  m and  n are  indices  locating  the  position 
of  a particular  subsection  over-which  an  inner  product  is  taken.  The  expres- 


' 


' 

i 
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slons  for  C(s) , ^mn's»  pmn's  anc*  q's  are  very  complicated.  They  have  been 
given  in  a previous  report  [22],  and  will  not  be  repeated  here. 

The  unknown  coefficient  matrices  {a(s)},  {B(s)>,  {-y(s)}  and  (J(s)}  can 
be  solved  from  Eq.  (32).  An  outline  of  the  method  of  solution  will  be  given 
in  Section  V.  Determination  of  these  coefficient  matrices  enables  the  calcu- 
lation of  f(r,s)  in  Eq.  (27)  from  Eq.  (31),  which,  in  turn,  leads  to  the  field 
Inside  the  cavity. 
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IV.  RECTANGULAR  CAVITY  WITH  SLOT  APERTURE 


PERPENDICULAR  POLARIZATION 


When  the  electric  field  of  the  incident  plane  wave  is  polarized  in 


a direction  normal  to  the  slot  shown  in  Fig.  1;  that  is,  if 


it  is  more  convenien't  to  use  an  H-formulation.  Instead  of  invoking  the 
induction  theorem,  we  apply  the  equivalence  theorem  [20],  [21].  The  H- 
field  discontinuity  at  the  aperture  in  Fig.  2(d)  is  supported  by  an  elec 


To  solve  this  problem,  we  start  with  the  following  Maxwell's  equations 


which  transform  to 


Let  A be  the  Laplace  transform  of  a magnetic  vector  potential  A such  that 


y eos  3x3y  x 


E ■ — r~—  A . 
z e s 9x3 z x 
o 


The  second-order  differential  equation  (40)  can  be  converted  into 
a set  of  first-order  equations  by  defining  new  quantities  0,  y,  and  O 
similar  to  those  in  Eqs.  (19)  to  (21)  with  F replaced  by  A . Instead  of 

X X 

Eqs.  (22)  to  (24),  we  now  have 


An  operator  equation  (29)  which  represents  Laplace  transformed  normalized 

state  equations  in  the  four  state  variables  A^,  u,  v,  and  w is  obtained 

when  operators  L and  P and  column  matrices  f(r,s)  and  e (s)  are  defined  as 

P 

in  Eqs.  (25)  to  (28)  respectively  with  F^  replaced  by  A^  and  replaced 
by  J^.  This  operator  equation  can  be  solved  by  the  method  of  moments  as 
before  by  choosing  suitable  expansion  functions  A . . , . . (r) , u,.  . (r) , 

v,.  . . . (r) , and  w,.  . ,.(r)  over  cavity  subsections.  We  can  write  the 


expanded  form  of  f(r,s)  inside  the  cavity  as  [22] 


f(r’S>  " ij.k  t5<1.J.k><S>  f(i.).k)<r)  + ®(i.J,k)<s)  f(i,j,k)<r) 


+ ,(l,J,k><9)  f(l,j.k)(r)  + 5'(i.3,k)(s)  <50) 

11  V w 

where  the  column  matrices  ^ ^(r),  ^ k)^’  ^(i  j k)  are 

the  same  as  those  defined  in  Eq.  (30), 


Ax(i,j,k)(r) 


f(i.3.W(r) 


and  a,  8,  f,  and  6 are  the  Laplace  transforms  of  the  expansion  coefficients. 

Analogously  to  the  parallel-polarization  case  treated  in  Section  III, 
we  can  expand  in  the  half-space  y > 0 as  a superposition  of  plane  waves: 


oo  oo 

f sV*-v 


-j (k  x+k  y+k  z) 
* y z 


dk  dk  , 
x z 


where 


jky  = /(s/c)*  + k£  + k£  (5 

and  the  new  quantity  g (,  ,k  ) can  be  determined  from  the  boundary  con- 

X z 

dltion  at  the  slot  [22],  There  is  a discontinuity  in  H at  the  slot  on 

z 

account  of  the  equivalent  electric  current  Jx  which  can  be  expanded  as 


V V(l.ny,k)(X>  °-  *’ 


4 
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A system  of  three  equations  are  obtained  by  matching  A^,  u,  and  at 
the  slot.  A combined  field  expression  which  holds  inside  the  cavity, 
in  the  y > 0 half-space,  as  well  as  in  the  slot  can  be  written.  This 
expression  is  extremely  complicated.  Suffice  it  to  say  that  when  it 
is  substituted  in  Eq.  (50)  and  inner  products  with  the  expansion  func- 

. 

tions  are  taken,  we  obtain  a matrix  equation  similar  to  Eq.  (32),  from 
which  the  unknown  coefficient  matrices  (a(s)},  (6(s)},  (y(s)}*  3°^  (S(s)} 
can  be  solved. 


V.  SOME  NUMERICAL  RESULTS 


In  this  section  we  outline  the  procedure  for  determining  the  unknown 

coefficient  matrices  (a(s)},  (8(s)},  {y(s)},  and  (6(s)}  from  Eq.  (32).  Al- 

, , , . , „ uv  „ vu  . W . vw  Fv  uv  vF  vu  w 

though  the  component  matrices  («.  , Z , Z , Z , p , p , p , p , p , 

mn  mn  inn  inn  mn  nm  mn  mn  inn 

and  p™)  representing  the  coupling  between  cavity  and  external  fields  appear 
mn 

highly  complex,  they  are  relatively  sparse.  Typically  an  equation  of  the 
following  form  is  obtained  from  Eq.  (32): 

[Z(s)  ] {a(s)>  = [H(s)  ] {-n  i \ (s)  } + (ic(s)  > , (55) 

s xQi,n^,K.; 

where  [Z(s)J  and  [H(s)]  are  square  matrices  containing  Z 's,  p 's,  and 

mn  mn 

q^'s  in  Eq.  (32),  and  (K(s) } is  in  general  not  the  same  as  (C(s)}.  From 
Eq.  (55)  we  have 


(&(s)}  = [Z(s)  ] 1 ( [H(s)  ] + <K(s)))  • (56) 

Similar  expressions  are  obtained  for  (8(s)},  {y(s)},  and  (6(s)}  . 

Let  s be  the  zeros  of  |z(s)|  or  the  roots  of  the  equation 
a 

det[Z(s) ] = 0 . (57) 

In  circuit-theory  terminology,  [Z(s)]  corresponds  to  the  system  impedance 
matrix  and  s^  are  the  natural  frequencies.  [Z(s)]  ^ can  be  expanded  in  a 
partial-fraction  form  as  follows: 


rz(s)]"1  = I j-~-  * (58) 

a a 

where  the  constant  square  matrix  [R^]  is  the  system  residue  matrix  at  the 

pole  s . [R  ] can  be  written  as  the  product  of  a natural  mode  vector  {Rm} 
a ot  a 

and  the  transpose  of  a coupling  vector  {R^}  [5],  [9]: 


[RJ  = (Rm}  (RC}T  , 

aJ  a a 


where  {Rm}  is  a solution  of  the  equation 


and  {R  } is  a solution  of 
a 


[Z(sa)]  {R“}  = 0 


[Z(sa)]T  {R®}  = 0 . 


A close  examination  of  the  composition  of  the  matrices  [Z(s)]  and 
[H<s) ] reveals  that  their  poles  coincide  and  therefore  cancel.  We  have, 
from  Eqs.  (56),  (58)  and  (59), 


{RmHRC}T 
{s(s)}  = l -1—2- 
s — s 

a a 


C[H(s)  ] {i*„(1>11  ,«(.))  + «(.)» 


Now  define 


[H(s)]  R,,  . v(s)}  + (K(s)}  = N(s)  (V  (s)}  , (62 

s x(i,n^,K;  o 

where  {V  (s) } is  the  excitation  vector  when  the  incident  wave  is  a pulse, 
o 

We  can  then  write  Eq.  (62)  as 

{Rn}{Rn}T  - 

(a(s)}  - l — — N(s)  {V  (s)} 

u s - s o 

a a 


where 


l Vs)  *(s) ' 


fi  (s)  = (Rc}T(V  (s)} 
a a o 


Li  * 

! i 


is  called  the  coupling  coefficient  [5].  We  note  that  N(s)  itself  may  con- 
tain poles  in  the  finite  plane,  but  this  fact  does  not  result  in  any  serious 


difficulty. 
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We  are  now  in  a position  to  write  the  expressions  for  the  field 


distributions  within  the  cavity,  ^rom  Eq.  (31), 


P*<W-S)  ‘ tjik  5(i,j,k)(8>  Fx<i,1,k)(*’y,*) 
■ ,S(S,)T 

which,  in  view  of  Eq.  (64),  becomes 


(66) 


Fx(x,y,z,s)  . J no(.)  {rJT  <•  - z,,)'1  N(«) 

= £ na(s)  vF(x,y,z)  (s  - s )-1  N<s>  . (67) 

a 

In  Eq.  (67), 

va(x,y,z)  = {Ra}  (68) 

is  a natural  mode  for  F . In  a similar  manner,  we  will  get 

E 

E (x,y,z,s)  = l n (s)  v y(x»y,z)  (s  - s )_J-  N(s)  (69) 

y a a.  a. 

and 

Ez  -1  - 

E (x,y,z,s)  = J fj  (s)  v (x,y,z)  (s  - s ) N(s),  (70) 

^ 01  Cl  Cl 

a 

E E 

y z — 

where  v J and  v are  the  natural  modes  for  E and  £ respectively, 
a a y z J 

For  numerical  computation,  a rectangular  cavity  of  dimensions 
4x6x2  was  chosen  and  the  slot  width  was  a/10  or  0.4.  Considerable  dif- 
ficulties were  experienced  in  the  determination  of  the  singularities 
(natural  frequencies)  sa.  After  having  compared  several  root-finding 
procedures  and  checked  meticulously  our  computer  programs,  we  obtained 
the  locations  of  the  singularities  in  the  upper  part  of  the  left  half- 


plane  for  the  parallel-polarization  case  with  E-formulation  as  shown  in 
Fig.  3.  Two  layers  of  singularities  were  found:  the  first  layer  was  very 
close  to  the  jw-axis  and  the  second  layer  had  a positive  slope.  The  singu- 
larities in  the  first  layer  fall  on  a smooth  locus,  and  they  appear  to  be 
irregularly  spaced.  No  method  exists  which  could  ascertain  the  effect  of 
numerical  noise,  subdivision  scheme,  etc.  It  is  known  that  numerical  results 
are  sensitive  to  round-off-errors  and  may  diverge  with  an  increasing  sampling 
density  [23].  Problems  of  numerical  anomaly  appear  in  even  the  simplest 
cases  [23],  and  precluded  the  determination  of  the  natural  frequencies  of 
some  well-defined  geometric  structures  such  as  oblate  spheroids  [24]. 

The  natural  modes  of  the  first  four  singularities  of  the  first  layer 

were  determined  and  combined  to  give  the  total  | E I at  two  locations  within 

z 

the  cavity.  Assuming  a step  excitation,  the  normalized  cavity  field  | [ /E^ 

is  plotted  in  Fig.  4 as  a function  of  time  at  the  point  (2,  -5,  1),  which  is 

(5/6) th  of  the  way  in  from  the  slot  aperture.  Figure  5 shows  a similar  plot 

at  the  point  (2, -3,1).  The  figures  indicate  that  the  maximum  value  of  [ E [ 

inside  the  cavity  is  in  the  order  of  one-thousandth  of  the  incident  field 

intensity,  E*,  for  parallel  polarization  (good  shielding  effectiveness), 
z 

Maxima  and  nulls  exist  as  t varies.  However,  a physical  interpretation  of 
their  spacings  is  difficult  because  of  the  complicated  nature  that  the 
penetrated  field  is  reflected  from  the  cavity  walls. 

Previous  studies  [25],  [26]  on  the  penetration  of  transient  electro- 
magnetic excitation  through  apertures  in  an  infinite  ground  screen  (without 
a cavity)  have  not  yielded  numerical  results.  Various  difficulties  were 
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encountered  in  attempts  to  obtain  self-consistent  pole  locations  because 
of  numerical  instability  [26].  The  problem  is  vastly  more  complicated 
when  there  is  a cavity  behind  the  aperture.  It  has  been  shown  [27]  that 
the  structure  of  a conducting  cylinder  within  a parallel-plate  region 
gives  rise  to  two  types  of  singularities;  namely,  poles  and  branch  points, 
in  the  complex-frequency  plane.  Whether  a cavity— backed  slot  has  branch 
points  or  multiple  poles  is  an  unanswered  question. 


[ 


| VI.  CONCLUSION 

In  this  report  the  transient  excitation  of  a cavity  through  an 
aperture  is  used  to  illustrate  the  state-space  solution  of  transient 
electromagnetic  problems.  Depending  on  whether  the  E-  or  the  H-formulation 
is  used,  the  effect  of  the  aperture  is  accounted  for  by  a magnetic  current 
in  invoking  the  induction  theorem  or  by  an  electric  current  in  invoking  the 
equivalence  theorem.  The  governing  second-order  differential  equation  is 
converted  into  a set  of  first-order  normalized  state  equations  by  defining 
three  new  state  variables  in  addition  to  the  appropriate  vector  potential. 
These  state  equations  subject  to  the  associated  boundary  conditions  are 
solved  by  the  method  of  moments.  The  cavity  region  is  first  divided  into 
subsections  and  the  fields  within  the  cavity  expressed  in  terms  of  appro- 
priate expansion  functions  with  time-dependent  coefficients.  The  fields 
in  the  half-space  outside  the  cavity  are  represented  as  superpositions  of 
plane  waves.  At  the  aperture,  the  cavitv  and  external  fields  are  properly 
matched.  Inner  products  are  taken  with  testing  functions  and  the  first- 
order  equations  are  converted  into  matrix  equations  containing  the  expansion 
coefficients  as  unknown  column  vectors.  Evaluation  of  these  expansion- 
coefficient  vectors  leads  to  the  determination  of  field  distributions.  The 
procedure  for  evaluating  a typical  coefficient  vector  by  the  singularity- 
expansion  method  has  been  outlined. 

This  solution  procedure  was  applied  to  the  problem  of  electromagnetic 
excitation  of  a rectangular  cavity  through  a slot  aperture.  Two  cases  were 
considered:  an  E-formulation  for  parallel  polarization  and  an  H-formulation 
for  perpendicular  polarization.  For  the  case  of  parallel  polarization,  some 

24 


significant  pole  singularities  (natural  frequencies)  in  the  upper  left 
complex  plane  were  formed  for  a sample  cavity.  The  corresponding  natural 
modes  were  determined  and  combined  to  give  the  electric  intensities  as 
functions  of  time  at  two  locations  in  the  cavity  for  a step-excitation. 
Both  the  shape  and  the  magnitude  of  the  penetrated  electric  field  appear 
reasonable.  However,  it  would  not  be  wise  to  claim  absolute  accuracy 
because  of  the  inherent  noise  and  instability  in  the  numerical  procedure 
and  because  of  a lack  of  a theoretical  proof  that  there  would  be  no  multi- 
ple poles  or  branch  points. 
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ture  can  be  accounted  for  by  a magnetic  current  in  invoking  the  induction 
theorem  or  by  an  electric  current  in  invoking  the  equivalence  theorem.  The 
governing  second-order  differential  equation  is  converted  into  a set  of  first- 
order  state  equations  by  defining  three  new  state  variables  in  addition  to  an 
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appropriate  vector  potential.  These  state  equations  are  solved  by  the  metho 
of  moments.  Two  cases  are  considered:  parallel  polarization  and  perpen- 
dicular polarization.  Some  significant  singularities  for  the  parallel- 
polarization  case  are  found  and  the  electric  intensities  as  functions  of 
time  at  two  locations  in  the  cavity  are  computed  for  a step  excitation. 
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